% CREATEMATRIXC
%
% deze functie maakt van een matrix A de juist testmatrix*.in en
% testmatrix*.out files aan voor het circulante geval

function creatematrixc(A,index)
delete(['../testmatrices/testmatrix' num2str(index) '.*'])

fileID = fopen(['../testmatrices/testmatrix' num2str(index) '.in'],'a');
fprintf(fileID,'c%2i \n',length(A));
for i=1:length(A)
fprintf(fileID,'%20.15e ',A(i));
fprintf(fileID,'\n');
end
fclose(fileID)
% out
fileID = fopen(['../testmatrices/testmatrix' num2str(index) '.out'],'a');
B = expm(circulant(A));
B = B(:,1);
fprintf(fileID,'c%2i \n',length(A));
for i=1:length(A)
fprintf(fileID,'%20.15e ',B(i));
fprintf(fileID,'\n');
end
fclose(fileID)